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By Bernoulli's law, an increase in the relative speed of a fluid around a body is accompanies 
by a decrease in the pressure. Therefore, a rotating body in a fluid stream experiences a force 
perpendicular to the motion of the fluid because of the unequal relative speed of the fluid across its 
surface. It is well known that light has a constant speed irrespective of the relative motion. Does 
a rotating body immersed in a stream of photons experience a Bernoulli-like force? We show that, 
indeed, a rotating dielectric cylinder experiences such a lateral force from an electromagnetic wave. 
In fact, the sign of the lateral force is the same as that of the fluid-mechanical analogue as long as 
the electric susceptibility is positive (e > eo), but for negative-susceptibility materials (e.g. metals) 
we show that the lateral force is in the opposite direction. Because these results are derived from 
a classical electromagnetic scattering problem, Mie-resonance enhancements that occur in other 
scattering phenomena also enhance the lateral force. 

PACS numbers: 



Photonic Bernoulli's Law? When considering a rotat- 
ing body in a fluid stream such as air, the body experi- 
ences a pressure gradient caused by the difference of the 
relative velocity of its motion to that of the fluid at var- 
ious points on its boundary. For example, an idealized 
tornado such as a spinning cylinder moves perpendicular 
to the streamlines of the fluid. The direction of motion is 
along the direction connecting the center of the cylinder 
to the point of maximum relative velocity. 

In a famous experiment, Michelson and Morley pQ 
showed that even if the earth were immersed in a fluid in 
motion, the speed of light would be constant relative to 
perpendicular directions. Later, the special theory of rel- 
ativity established the constancy of the speed of light re- 
gardless of observer's relative motion to the light source. 
Here, we ask to what extent can a stream of photons re- 
sembles a stream of massive fluids? In particular, if one 
considers a stream of photons (classically described by 
Maxwell's equations) as a fluid in motion and places a 
rotating dielectric body in it, one might naively expect 
that no Bernoulli-type force would be experienced by the 
body since the relative speed of light is the same on both 
sides. Here we show that such a force is experienced by 
the rotating body, though the cause is the asymmetry of 
the scattered field \§\ from the dielectric, by which a net 
force is imparted to the rotating body. 

Cylindrically Rotating Dielectric- The exact electro- 
magnetic constitutive equations in a medium moving at 
velocity v, discovered by Minkowski [3J, are 

D + vxH/c = e(E + vxB/c) (1) 
B + Exv/c = //(H + Dxv/c), (2) 

where E, D, B and H are the usual electromagnetic 
fields, c is the speed of light in vacuum and e is the electric 




Figure 1: Light scattering from a rotating dielectric cylinder. 



permittivity in the rest frame and [i is the magnetic per- 
meability in the rest frame. These equations presuppose 
uniform motion of the dielectric, where special relativ- 
ity is sufficient. For accelerated dielectrics, the equations 
become more complicated; however, for rotating bodies 
with axial symmetry, the body in motion has the same 
shape as the one in the rest frame and it has been shown 
that the same equations would apply [1HB]. This asser- 
tion has been successfully used in applications [B] and 
was later proved rigorously by Ridgely [5], who showed 
that the general relativistic treatment for uniformly ro- 
tating dielectrics with axial-symmetry, to first order in 
v/c, gives Minkowski's results (Eqs. [I]and[2]). 

In the limit where v/c is small, Tai considered the scat- 
tered field of a plane wave incident upon a uniformly ro- 
tating dielectric cylinder with angular speed f2 [6]. We 
begin by reviewing Tai's derivation of the scattered field 
and then we use these fields to compute the force. As 
depicted in Fig. [T] the velocity of the rotating body is 
v = rtrcj) at a radius r, the radius of the cylinder is 
denoted by a, and the E field of the incident wave is 
assumed to be polarized in the direction of the axis of 
the cylinder (which we take to be z). Derivations of key 
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equations are provided in the appendix. 

We solve the scattering problem by standard technique 
of expanding the field in each region in basis of Bessel 
functions J n and then matching boundary conditions at 
the interface. In this basis, an incident z-polarized plane 
wave propagating in the +x direction with amplitude Eq 
(see Fig. [IJ is given in polar (r, <f) coordinates by 



Ei = E exp {ikgr cos (/>} = Eq i" J n (kor) ex P (in^>) , 

n— — oo 

(3) 

where Eq is the amplitude, kg = lo/c is the wave num- 
ber in vacuum, u is the frequency in the time-harmonic 
oscillating field e~ lujt . The scattered and "transmitted" 
(interior) fields, respectively, can be written (using the 
Hankel function H n x ^ = J n + iY n ) with to-be-determined 
coefficients a n and 8 n : 



+00 

E s = E J2 a n i n H^(k r)exp(in^). (4) 

71= — OO 
+OO 

E t = E Q Pni n Jn {j n r) exp (incf>) , (5) 



where j n is defined by 



K = Ho(e- eo) 

7- 



(6) 



2 - k 2 - 2nwA' = /c 2 I 1 I . (7) 



The total field is therefore E = E z z = (Ei + E s )z and 
the magnetic field in the vacuum regions is given by H = 
-J— V x E. 



The unknown coefficients /?„ and a n are found by re- 
quiring continuity of E z and iJ^ at r — a, yielding 



J n (koa) + a n H^ (k a) = P n J n (j n a) (8) 



fc 



J' n (k a) + a n H^ (k a)' = /3 n j n J' n (j n a) , (9) 



where the prime on J and if W denotes derivative with 
respect to the entire argument. Solving for a n gives 



Jn (Po) [Jn-l (Pn) ~ Jn+l (/Jn)] - V~Jn (Pn) [Jfi-1 (Po) - Jn+l (po)] 



M 1} (PO) [Jn-l (Pn) - Jn+l (Pn)] ~ ^ Jn (Pn) (p ) - (p ) 



(10) 



where p = fc ° an d Pn = 7n a - For 51 ^ 0, the rotation 
breaks the y — mirror symmetry leading to asymmet- 
rical scattering a n 7^ a_„ as shown by Tai [6J. If Q = 0, 
then a„ = a_ n and Eq. [4] reduces to symmetrical scat- 
tering. 

Force imparted to the Rotating Cylinder.- The asym- 
metry in the momentum transport by the scattered field 
should manifest itself as a lateral force on the dielectric. 
This force can be computed by integrating the Maxwell 
stress tensor over a closed surface around the object. Be- 
cause we only evaluate the stress tensor in vacuum, we 
avoid the well-known difficulties that arise in defining the 
stress tensor inside the material [2], nor does the rotation 
affect the vacuum stress tensor. The stress tensor in SI 
units is 



++ 



W = e E ® E + /i H ® H 

- (e E 2 + fi Q H 2 ) (x<g)x + y®y + z<g)z 



(11) 



where the hatted quantities are unit vectors. To calculate 



the force on the cylinder in any direction ri on the plane 
at a fixed radius tq, we evaluate 



2^ Jo 



dt (h r dr d(p 8 (r — r ) {ri • xr ■ r} 



(12) 

where f = cos </>x + sin </>y and the time average is taken 
over a full period to obtain a real force. For our polar- 
ization, E x = E y = H z = 0. The force in y direction 
is 



Fy = r 



e |^| 2 + p |H| 2 



sin</> (13) 



W D „,u, J. , Po \H y \ 



+ =iRe cosy 2 



sm ( 



Using orthogonality conditions on e 1 ™* Eq. 13 can be 
integrated analytically. In the figures below we use the 
dimensionless force where P = 2a \Eq\ 2 is the 
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Normalized Force vs. Rotational Frequency for*. = 0.01 



Normalized Force vs. Rotational Frequency for X = 0.01 
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Figure 2: Force vs. rotational frequency for e/en = 10 and 
P = 2a\E \ 2 .f^. 
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incident power on the scatterer's geometric cross section. 
This is a convenient normalization because, for light in- 
cident on a perfectly absorbing flat surface the force is 
exactly —, so this normalization gives a measure of the 
lateral force relative to the incident photon pressure. Fur- 
thermore, we use the dimensionless angular frequency 
^p, which is the ratio of speed of the cylinder bound- 
ary to that of light. 

Fig. |2]plots the force vs. angular frequency for e/eo = 
10. As discussed below, whenever e > eo, corresponding 
to a positive electric susceptibility = c/cq — 1, there is 
a force of sign analogous to Bernoulli's law, where > 
(counterclockwise rotation) gives a force in the positive 
y direction. Clearly, ft < gives the same magnitude 
of the force in the opposite direction as expected from 
the symmetry of the problem and in accordance with the 
fluid-mechanical analogy. 

In the case of a perfect conductor, j3 n is zero (see Eq. 
|5| and therefore Eqs. [8] and [9] do not give an asymme- 
try with respect to n for a n . Consequently, in this limit 
there is no lateral force. Intuitively, because a perfect 
conductor allows no penetration of the electromagnetic 
fields, the fields cannot "notice" that it is rotating or be 
"dragged" by the moving matter. However, for imperfect 
metals (finite e < 0) there is some penetration of the ra- 
diation into the material which results in a lateral force. 
Interestingly, in the case of e < 0, and in fact whenever 
e < eo (negative susceptibility), the force is in the oppo- 
site direction of the force for e > eo dielectrics (see Fig. 
|3|. The reason is an immediate consequence of Eq. [6] 
For e < e and ft > 0, K becomes negative and the phe- 
nomenology, looking at 7„ in Eq. [7] become equivalent 
to the case of e > eo and f2 < 0. The same relationship 
between the sign of the force and the sign of Re e — eo 
holds for complex e as long as |Ime| <C |Ree|, whereas 
for large | Ime| we observe a similar relationship with the 
sign of the imaginary part. 

Lastly, we investigate the dependence of the normal- 
ized force on 2nuja/c = a/Ao, varying the vacuum wave- 
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Figure 3: Force vs. Rotational frequency for various e/eo and 
P = 2a\E \ 
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Normalized Force vs. J. for (a D)/c=3.34e-04 and a=0.01 
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Figure 4: Mie resonances: Force vs. An for a — 0.01 



length Ao (see Fig. |4J). For Ao -C a the scattering ap- 
proaches a ray-optics limit, while for Ao 3> a it is in the 
Rayleigh-scattering (dipole approximation) regime [S]. 
For Ao ~ a, the force spectrum becomes more interesting 
due to the presence of Mie resonances [S] ■ 

Discussion and Future Work.- The forces obtained 
here are only a fraction of the incident radiation pressure 
and seem to require infeasible rotation rates, but we ex- 
pect that they can be resonantly enhanced by techniques 
similar to those that have been used by other authors 
to enhance scattered power for a given particle diame- 
ter. Mie resonances are already visible in Fig. [4] but 
much stronger resonant phenomena can be designed by 
using multilayer spheres that trap light using Bragg mir- 
rors and/or specially designed surface plasmons, and one 
can even obtain "superscattering" by aligning multiple 
resonances at the same frequency jTU]. Such enhance- 
ment mechanism, in combination with recent progress in 
generating rotating particles (of graphene) at near-GHz 



4 



[TT], may permit the future experimental observation 
and exploitation of optical "Bernoulli" forces. 
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I. APPENDIX 

First let a, /i, e be the conductivity, permeability and 
permittivity respectively and define A = (e/i — eo/Lto)v, 
where subscript corresponds to quantities in the vac- 
uum. The electric and magnetic field vectors satisfy 



V x E = - — (uH — A x E) 

dt v ' 



(14) 



V x H = a (E + /iv x H) + — (<jE + A x H) .(15) 

at 



To determine the proper expression for the transmitted 
wave, we use harmonically oscillating fields to reduce Eqs. 
to (neglecting O ^(f/c) 2 ^ terms) 



14 



and 
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(V + iwA)xE = iufiH. (16) 
(V + iu)A - cr/iv) x H = (a- iuje) E (17) 

Looking at Figure [T] the incident wave is given by 



Ei = E exp {ikor cos </>} = E i n J n (kor) exp (in(/>) , 

n— — - oo 

(18) 

where fco — w / c is the wave number in vacuo; u being 
the frequency in the time harmonic oscillating field e lut . 
The scattering field can be written in the form 



rnc 



(19) 



n— — oo 



To solve for the transmitted field inside the dielectric we 



subject Eqs. 16 and 17 to the particular form of the 
velocity which is independent of z. The parameter A is 
then a function of r alone, 



A =^ (e - eo) v = Krcj) 



(20) 



where K = fioe (l 



= 



mil 



with m — 1 



and 



20 



into Eqs. 



16 



and 17 and 



c = 1//j,q€. Substituting Eq. 
eliminating H we obtain a differential equation for E z 
which is the only component of the electric field inside 
the dielectric cylinder. 



ld_ 

r dr 



,dEz 
dr 



1 d 2 E z 



90 2 



-2iwK- 



.dE z 



-k 2 E z +0 



= 0, 



(21) 



He) 



where k 2 = 
tions for E 

understanding that the results are accurate to first order 



To solve let us seek separable solu- 
F (r) e"^ and below we drop O ((-) 2 ) by 
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The function F (r) then satisfies 



1 d ( dF 



2tiljK — k } F = 



r dr V dr 



If we introduce 7^ — k 2 — 2nuK = k 2 (l 



2nmf2 > 



then the proper set of radial functions to describe the 
field inside the rotating cylinder is 



fc 



J n (k a) + a n H^ (koa) = j3 n J n (j n a) (23) 
J' n {k Q a)+a n H^ (k a)'} = /3 n7n J' n ( 7n o) ,(24) 



where the prime on J and -f/^ 1 -* denotes derivative with 
respect to the entire argument of these functions. The 
solutions for a n and /3 n are 



F(r) = J n (7„r), u = 0,±1,±2, • • • . 

The complete expression for the transmitted field can be 
written in the form 



+00 



E t = E ^2 Pni n j n {l n r) exp (i 



(22) 



n— — 00 



By matching E z as defined by Eqs. [18j[T9]and[22j and the 
component of the magnetic field at the boundary r — 
a one obtains the following two simultaneous equations: 



Jn (po) J'n On) - ^ J n (Pn) J' n (po) 

(25) 



H ( n ] (Po) J'n (Pn) ~ ^Jn (Pn) H™' ( Po ) 

k 



Jn (p ) Hn ] ' {po) - J'n (Po) ( P o) 



(-26) 



7n H n 1] (po) J' n ( P n) ~ ^Jn (p n ) H^ h ( P o) 

where p Q = k a and p n = j n a. Using identities J' n (x) = 
\ J n _i (x) - J n+1 (x)], J' Q (x) = -Ji (x) and H^' (x) = 
\ [^1, (x) - Hi% (x)] , H^' (x) = -H^ (x), f3 n can 
be eliminated to give 



Jn (po) [Jn-l (pn) - Jn+l (Pn)] ~ ^J n (p n ) [Jfl-1 (Po) - Jn+1 (Po) 



a n = 



Hn 1} (Po) [Jn-l (Pn) ~ J n+ 1 (Pn)] - ^ Jn (Pn) (p Q ) - ( Po ) 



(27) 



The numerical value of a n 7^ a- n for f2 7^ because 
in this case p n 7^ p_„ and j n 7^ 7_ n hence the scattering 
field has an asymmetrical part with respect to the direc- 
tion of incidence, = 0. When ft = ==> a n = a_„ 
and Eq. [19] reduces to the well known results. 

The total field therefore is E = E z z = (Ei + E s )z. 
Below we suppress kor as the argument of the Bessel 
functions unless stated otherwise. 



Let rib =cos')/ , x + sin^y be any unit vector, then Eq. 
[PT] evaluated at the radius ro, reads 



*W 



r <f> d<j> 
Po 



Po 



\H\ 



cos (ip - </>) (30) 



Re (H*H y ) sin (ijj 



+ OO 



E 



(28) 



n— — 00 



Further in free space we have H 



^— V x E which 

ILJUq 



gives 



Po \Hx\ 



cos 4> cos ip 



Po \H y \ 



sin <fi sin ip 



In particular we are interested in ip — ^ to calculate 
the transverse force 



H 



— V x Ez= — 



iujpo 
1 fldE z 



iujpo \r o<p 
1 ( 1 dE z 
iujpo 



cos 



sin < 



lujpo \r 
, dEz 
Or 
dE* 
dr 



ldE z , dE z 

r- -^—<j) (29) 



dr 



sine 



COS( 



e \E\ 2 +p \H\ 2 



+ ^Re (#;#,) COS 



where _E = (^ + _E S ) is given by Eqs. [3]and[4j 



(31) 
(32) 



II. EVALUATING THE INTEGRAL 



where 



Here we evaluate \E\ , \H\ , \H y \ and Re(H*H y ) as 
they are useful for calculating the force below (Eqs. [30| 
and 13). The key equations are 



H 



E Q i n {Jn + a n H^} e m « 



n— — oo 



H x ± + H y y 
1 (ldE z 



iuj^o \r d(f> 
1 (ldE z 



, dE z . 
cos <p H — - — sin0 | x 
or 



dE z 
dr 



cos y 



H„ = 



E 



+ 00 



lUJLLn *~ 
r n——oo 



1 1 



in<f)-n 



E 



4WUn '- 

r n=—oo 



' 1 



I 4- n> 



I 4- (T W 



COS< 



Sill — 



fc 

2 

Ay, 



--ff. 



Ci) 



n+1 



Jn—1 Jn+l "I" C^n 



COS « 



as well as the orthogonality relations 



We proceed 



6e i(n-m)0 gin( 



be i(n-m)4> gin 3 ( 



ie i(n-m)^ gin 2 ^ cog( 



^ ( "- m) *sin (icos 2 . 



6e i(«-™)0 cos 3, 



(-in) (<5„ :m _i - <5„, m+ i) 

(^n,m-3 — ^n,m+3 — 3(V m _i + 35 n ^ m +i) 
(^n,m— 3 ~t~ ^n,m+3 ^n,m— 1 ^n,m+l) 
(^n,m— 3 ^n.m+3 "T" ^n,m — 1 ^n,m+l) 
(4) ^™> m ~ 3 ^™> m +3 + 35„ im _i + 35„, m +i) . 



?7T 

T 

— 7T 

— in 



/+00 „ 
#|£ 2 | 2 sm</> = S 2 £ *d^ i ( B - w J*sin^(-ir^{j n + anH( 1 >}{j ro + a^fl-W} 
n,m— — 00 
+00 

= £ 2 ^ (-I)™ I n+m [j n + a n ff«} { J m + C4i4? } } HO {<J„, m -l - S n , m+1 } 



n.m— — 00 



Further § d<fi \H X \ 2 sin <f) — £ dcf) (H*H X ) sin <f>, similarly 



for 



H v I sin ( 



(j) d(j) |ffa;| 2 sin( 



\U}fJ, 



2 +00 „ 
^ {-l) m i n+m f, 



> sin me 



i(n—m)4> 



n.m— — 00 



\j +a* HV> 

•Jm 1 {JL m 11 m 



COS(j> + 



k 



t _ 7 Mr,* - \ 



sin i 



in 

■ I — 

r 



COS + — 



J n -1- Jn+l+On (H^-H^) 



sine 



which using the orthogonality relations yields 



I dci>\H x \ 2 sin <f> = (^-\ ]f (-1) 



nm 



(^n,m— 3 ^n,m+3 H~ <^n,m— 1 <^n,m+l) 



+ 



Jm-l - ^m+1 + "m (#m-l _ #m+l) ^"-1 _ + a ™ (#n-l - ^n+l) 



' "4" ) (^«>™-3 — <5n,m+3 — 3(5„ !m _i + 35„ jTO +l) 



+ 



27 



-m (j m + cf m Hg>) - J n+1 + a n - itf^)) 

+n (j„ + a„ff«) (j ro _x - J m+1 + < (i^ - H^ + i)) 



I ) (#n,m-3 + $n,m+3 — <^n,m-l — <^n,m+l) 



Similarly 



& d4>\H y \* 4*4, = (-^Y £ (-i) m i»W 



) sin me 



i{n—m) 



n,m—— 00 



sin 

2 



(7 n -|- cx n H^ l ^ 



sine 



7 7 +rv* (V (2) - 77 (2) ^ 

Jn-1 - Jn+1 + OL n (h^Ix - 



COS( 



COS ( 



which using the orthogonality relations yields 
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f d(f> \Hy 



2 . 
Sill ( 



/ p \ 2 + °° 

(£) £ (-')™'" + ™ 



n.m— — oo 



mn 



,ro— 3 — ^n,m+3 — 3(5 rl;m — i + 3<5 n . m +i) 



+ (jm-1 - Jm+1 + a* m (jl^-l ~ H m+l)) (jn-1 ~ Jn+l + «n (#1-1 - #l+l)) 



' I ~ ^~ ) (<$n : m-3 — <5n,m+3 + $n,m-l — ^n,ro+l) 



+ 



ifco 
27 



(i) 
+i 



m (j ro + a* m H^ (j„_! - J„ +1 + a„ (fl^ - fl£ 
-n (j n + (j m _x - J m+1 + < (tf^ - ff^)) 



I ~y J (^n,m-3 + ^n,m+3 — <^n,m-l — ^n : m+l) • 



Lastly we need § d(f>Re (H*H y ) cos < 
i <f d<j)H*H y cos + c.c. 



jdc»H* x H y 



COS( 



2 +oo „ 

53 {-\) m i n+m L 



> cos 0e 



i(n—m)4> 



n,m— — oo 



in 
I — 

r 



— [ Jm + a* m H ™ COS + y J m _i - J m+ i + a m (h^-i - H , 
J n + a n H^] sin<£ - ^ [j„_x - J n+1 + a n (h^ - JT#i) 



(2) 
m+1 



which using the orthogonality relations yields 



\m -n+m 



/( E \ 
d<f>H*H y cos 4> = 
V MO/ 7), j 777, — — OO 

(j m + (■/„ + an^ 1 ') - f (Jm-i - J m +i + a* m - R 



mn 



(2) 
m+1 



X (yJn-l — Jn+1 + 0> n (^H^\ - H^_^ ^— — ^ (<5„, m _ 3 — <5 n , TO +3 + <5„ !m -l - 5 n ,m+l) 
+ ^ (J m + O^) (Jn-1 Jn+1 + a n (H^ H^)) 

X (4) ^™' m ^ 3 ^ n < m + 3 + 3<5„, m _i + 3<5„, m +i) 



■ I -y ) (Sn,m-3 + <5n,m+3 — <5n,m-l — <5n : m+l) • 
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All of the above integrals were checked against numer- 
ics before calculating the cumulative effect that appears 
in the force Eq. [13] The total force was also checked 



against numerical experiments. In all cases agreements 
were found with errors of order O (lO -26 ). 



